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Abstract 

The quasi-static solutions of the matter density perturbation in F(R) gravity models have been 
investigated in numerous papers. However, the oscillating solutions in F(R) gravity models have 
not been investigated enough so far. In this paper, the oscillating solutions are also examined 
by using appropriate approximations. And the behaviors of the matter density perturbation in 
F(R) gravity models with singular evolutions of the physical parameters are shortly investigated 
as applications of the approximated calculations. 
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I. INTRODUCTION 



In order to explain the observational results about current accelerated expansion of the 
Universe, we need to modify the Einstein equation or introduce the anisotropic metric. The 
standard way to explain the accelerated expansion of the Universe is by introducing the new 
energy called "dark energy," or revising the left-hand side of the Einstein equation to ob- 
tain "modified gravity." The most popular model of these is the Lambda cold dark matter 
(ACDM) model because of its simplicity. The ACDM model is given by introducing the 
contributions from the cosmological constant and the cold dark matter into the Einstein 
equation. However, an extraordinary fine-tuning for the value of the cosmological constant 
is necessary to explain the observational results, because the typical scale of gravity Mp\ 
is much larger than the observed scale of the critical density ~ 10~ 3 eV. Therefore, there 
have been suggested many dark energy models and modified gravity models to relax the 
fine-tuning problem. As an example of modified gravity models, there are F(R) gravity 
models[l-5]. F(R) gravity models are described as an extension of the Einstein-Hilbert 
action to have higher power terms of the scalar curvature R. Today, viable F(R) grav- 
ity models cannot be distinguished from the ACDM model by the observational data of the 
background evolution, so it is necessary to evaluate carefully the perturbations in these mod- 
els. The perturbations in F(R) gravity models had been also examined before the discovery 
of the current accelerated expansion of the Universe . However, we need to reconsider 
the perturbations because the matter can be negligible during the inflationary era of the 
early universe, but the matter is an essential component of the late-time Universe. And we 
are interested in the evolution of the matter density perturbation. In this paper, not only 
quasi-static solutions, but also oscillating solutions of the matter density perturbation are 
investigated by analytic calculations. By numerical calculations, this oscillating behavior 



of the matter density perturbation has been found in L. Pogosian and A. Silvestri|8| and 
S. Carloni, E. Elizalde, and S. Odintsov 9)]. The behavior of the matter density perturba- 
tion in F(R) gravity models with singular evolutions of the physical parameters are also 
investigated. 

In Sec. II, the matter density perturbation is investigated by using some approximations 
and as applications, the perturbation is considered in the models which have singular behav- 
ior of some parameters. Section III begins from the derivation of the differential equation 
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without using any approximations 10|, and the behavior of the matter density perturbation 
is investigated while using the case analysis. The popular viable models of F(R) grav- 
ity are considered in Sec. IV. And concluding remarks are given in Sec. V. The units of 
&b = c = h — 1 are used and gravitational constant 8irG is denoted by k 2 = 8%/Mp 2 with 
the Planck mass of M P1 = G~ x l 2 = 1.2 x 10 19 GeV in this paper. 



II. PERTURB ATIVE EQUATIONS AND APPROXIMATED CALCULATIONS 
A. Background equations and linearized equations 

Here, we adopt the following action as a form of F(R) gravity models: 

S = 2^ / d * x ^ \- R + K R )] + 5matter ' ^ 

where / is an arbitrary function of the scalar curvature R, and f(R) represents the deviation 
from the Einstein gravity. When we use the spatially flat Friedmann-Lemaitre-Robertson- 
Walker metric, ds 2 = a 2 (r])(dr] 2 — J^ i=1 dx l dx l ), the Friedmann-Lemaitre equations are writ- 
ten by 

^(1 + f R )-\{R + f)- ™f R = (2-2) 
- 2 (H' + 2H 2 )(l + f R ) - l(R + f)- - 2 (Hf' R + f^) = k 2 w P , (2.3) 

where R = 6a~ 2 (7i' + H 2 ), f R = df(R)/dR, and the prime represents the differentiation 
with respect to conformal time rj. p is the energy density of the matter coming from the 
variation of Setter an d w is the equation of state parameter expressed by w = p/p. The 
Hubble rate with respect to conformal time % is defined by T-L = a' /a. The equation of 
continuity is written as follows: 

pi + 3(1 + w)Up = 0. (2.4) 

The above equations fl2.2l) - fl2.4p are background equations of the Universe. 

On the other hand, (00), (ii), (Oi) = (zO), and (ij), i ^ j elements of the linearized 
Einstein equation in the Fourier space are represented as follows, respectively, when we use 
the metric of the Newtonian gauge, ds 2 = a 2 (7/)[(l + 2^)dr] 2 — (1 + 2\1>) J2i=i dx l dx l ]\ 

(1 + f R ){ - A; 2 ($ + - 3H(& + + (3H' - 6^ 2 )$ - 3H'^} 
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+ f' R {-m$ + 3ft$ - W) = K 2 pa 2 S, (2.5) 

(i + f R ){<s>"+y" + m{® + + m'§ + in' + 2u 2 )^} 

+f' R (3H$ - ft^ + 3$') + /£(3$ - W) = c 2 K 2 pa 2 5, (2.6) 
(1 + + W + ft($ + *)} + /^(2$ - *) = -K 2 pa 2 {\ + (2.7) 

$ - * - ^7t^7tt{3^" + 6(ft' + ft 2 )$ + 3ft($' + 3tf') - fc 2 ($ - 2^)} = 0. (2.8) 
a 2 (l + Jr) 

Here, we define the sound speed by c 2 = 5p/5p and the matter density perturbation by 
5 = Sp/p, to express Eqs. (j2.5|) - (l2.8p . From the perturbation of the equation of continuity, 
V^T^ = 0, we have 

3tf'(l +w)-6' + 3H(w - c 2 s )S + k 2 (l + w)v = (2.9) 

and 

c 2 

$ + — — <5 + w' + ftv(l -3w) = 0. (2.10) 
l + w 

We treat the equation of state parameter and the sound speed of the matter as w = c s = 
by assuming the era after matter dominance. 



B. Sub-horizon approximation 

In the following, we derive the differential equation of the matter density perturbation 
by using sub-horizon approximation. Oscillating solutions cannot be, however, derived by 



such a procedure as shown in the fc-essence model [ll|, so we will discuss them in the next 
subsection. Formally, sub-horizon approximation is given as d/drj ~ ft <C k. It contains two 
approximations, which are small scale approximation ft <C k and quasi-static approximation 
d/dr] ~ ft. From now on, we apply sub-horizon approximation to Eqs. (I2.5p - (l2.8p to deduce 
the equation of the matter density perturbation; then it is necessary to be careful of the 
order of calculations. Unlike in the case of the exact calculations, the procedure to derive the 
equation of the matter density perturbation is not unique in the approximated calculations. 
For example, we can choose at least two procedures to derive the equation with respect to 
S when \fKR.\k 2 /a 2 <C 1. The first procedure is as follows: 

1. Combining Eqs. (I2.9p and (I2.10P to erase v. 

2. From Eq. (I2.8p . reading a relation $ « ty. 



3. By using a sub-horizon approximation, we obtain — 2k 2 (l + ~ n 2 pa 2 5 from 
Eq. (|23P. 

4. Gathering the above equations gives the equation only described by the parameter 5. 
The other procedure is as follows: 

1. Combining Eqs. (12. 7p and ( 12. 9 p to erase v. 

2. From Eq. d 2 . 8 j) . reading a relation $ « ^. 

3. By using a sub-horizon approximation, we obtain — 2k 2 (l + ~ K 2 pa 2 5 from 
Eq. (T23D. 

4. Gathering the above equations gives the equation only described by the parameter 5. 

The only difference between the two procedures is the first step. Because of it, the derived 
equations, however, become quite different. So we should be sure of which is the correct 
equation. The answer is obvious. The approximated calculations we considered by now are 
only correct at the leading order, but the cancellation of the leading order happened in the 
second procedure, so that the correct equation is the first one: 

P+W -WT73' = °- (2 - n) 

Similarly, when \fKR\k 2 /a 2 ^> 1, we have the following equation by using sub-horizon ap- 
proximation: 

5" + H8'-£^8 = 0, (2.12) 
(4 + jr) 

because a relation $ « 2\1/ can be read from Eq. (12.8j) . 



C. The approximation of ultra-relativistic variation 

As referred in the last subsection, oscillating solutions cannot be derived by sub-horizon 
approximation, so we use the approximation of ultra-relativistic variation d/dr] ~ k and 
small scale approximation T-L k instead. To obtain the equation of S, we first assume a 
relation $ w — ^ from the leading terms of Eq. ( 12. 6p . And then, Eq. (12. 8p gives 



If \ fRR\k 2 /a 2 <C 1, there is an obvious solution so that the F(R) gravity models which satisfy 
the condition \fRR\k 2 /a 2 <C 1 do not have oscillating solutions of \& and 5 within the form 
of e lhv . On the other hand, if \fnR\k 2 /a 2 3> 1, \& should satisfy the equation, ty" + k 2 ^ « 0. 
The solutions of it is, 

ijf = C ie - ife " + C 2 e ifo? , (2.14) 

where C\ and C<i are integration constants. Then, combining Eqs. (12. 9 j) and (I2.10p gives a 
relation between 5 and \& as 5" ~ 2\1/", so we have solutions of 5: 

6 = 2C^- lkr] + 2C 2 e ife " + C 3 r/ + C 4 , (2.15) 

where C3 and C4 are integration constants. However, the behavior of the oscillating solu- 
tions are not only determined by the leading order terms. Therefore, calculations without 
approximations are necessary for correct evaluations. 



D. The behavior near the singularity 



It is known that there are several future singularity scenarios of the Universe 12j, |13] 



• Type I ("Big Rip") : For t — > t s , a — > 00, p c ff — > 00, and \p c s\ 00. This also includes 
the case of p e g, p e s being finite at t s . 

• Type II ("Sudden") : For t — >■ t s , a — > a s , p c s — > p s , and \p c s\ 00. 

• Type III : For t — > t s , a — > a s , p c s — > 00, and \p e s\ — > 00. 

• Type IV : For t — > t s , a — > a s , p c s — > 0, |p e fr| — > 0, and higher derivatives of H diverge. 
This also includes the case in which p e g (p e s) or both of p e s and p e ff tend to some 
finite values, while higher derivatives of H diverge. 

And lately the little ripjl^] scenario is discussed, where any finite singularity does not appear 
but causes the destruction of the bound system. In the following, we will investigate the 
behavior of the quasi-static solutions of the matter density perturbation near these singu- 
larities. However, we do not consider Type I, Type III, and the little rip scenarios, because 
the approximations in these scenarios are equivalent to the superhorizon approximation. 
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Type II 



In the case of Type II, the Hubble rate H = a(t)/a(t) is finite, but the derivative of the 
Hubble rate H is singular, so the approximations k 2 , U 2 -C \U'\ are applied in the following. 
And the quasi-static approximation d n ~ U will be used. We first consider Eq. (j2.8p to 
obtain the equation with respect to 5, because the equation ty" « — 2"H'$ can be read by 
considering relations W ~ (U^>)' ~ U'^f. And 2$ w \I/ is assumed from Eqs. (12. 7p and 
( 12. 9p . Then, the equation of ty" + W^f = is obtained. The solutions of the equation 
can be written as \& = const. /a + const.e~ a , then ty' = —U^>. By using Eq. f 12 . 5 j) and 

= —U^f, we have the following equation: 

*» » 2tf (2.16) 

Here, the assumption 2$ « \l/ is justified because Eq. (12.161) implies that \1/ 3> |K 4 5/(fc 2 'H / )|. 
To erase the parameter v from Eqs. (12.91) and (12.101) gives the following equation: 

5" + US' - W - %WH' + fc 2 $ = 0, (2.17) 

so that the substitution of Eq. (I2.16P into Eq. (I2.17P gives the equation of the matter density 
perturbation, 

Equation fl 2 . 1 8 [) indicates that the behavior of the matter density perturbation changes near 
the singularity because Eq. (I2.18P is different from Eqs. (12. lip and (I2.12p . 



Type IV 

In case of type IV, the approximations k 3 , U 3 , U\U'\ <C \U"\ should be applied, and 
from Eqs. (Q and Q U 2 \f RR \, U A \f RRR \ < 1 are satisfied. Then, Eq. flZSJ gives $ » ^. 
Equations (1231) . fl2~7]) . and (Q give 

(2.19) 

The substitution of Eq. (12 . 19[) into Eq. (I2.17P will give the third order differential equation 
with respect to 5, but it is not explicitly expressed here. 
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There is a notable literature written by A. de la Cruz-Dombriz et al. |10|, which expresses 
a strict derivation of the differential equation of the matter density perturbation in F(R) 
gravity models. However, it is too difficult to analyze the time evolution of the matter 
density perturbation exactly because the differential equation is so complicated. We use 
some approximations in the following. We then need to be careful to the order of applying 
approximations, because their order is not commutable. 

III. ANALYSES STARTED FROM THE FOURTH ORDER EQUATION 
Case I 

First, the differential equation of the matter density perturbation Sp/p when we use 
small scale approximation H 2 <C k 2 is derived by using Eq. (31) and the Appendix of [k]] 
as follows: 

6"" + « (3 + f * + 0(H 2 /t 2 )) 6>" + « ! + 0(1)) 6" 

+ tf g + 0(1) ) i -_^( 5 ^_ + 0(1) ) 4 = „, (3,) 

where we have neglected sub- leading terms under the approximation T-L 2 -C k 2 . If we assume 
quasi-static solutions of 5, then we obtain an equation from the terms proportional to k 2 in 
Eq. (JO): 

d 2 5 (I 3 \ d8 2 . n . . 

n + Wes ^m<$ = 0. (3.2) 

dN 2 \2 2 ctt J dN 1 + f R m V ; 

Here, we have defined w eS = -2H/(3H 2 ) - 1, Q m = n 2 p/(3H 2 ), and N = In a. And H 

means the standard Hubble rate H = a(t)/a(t). In the matter dominant era, Eq. (13.21) 

becomes 

" liS 2 « = 0. (3.3) 



dN 2 2dN l + f R 

The growth factor of 5, f = <iln \5\/dN, is given by / = — 1/4 ± v^33/4 when \ f R \ <^ 1 from 
Eq. (13. 3p . Therefore, the growth factor of this case in the matter dominant era is bigger than 
that of ACDM model, / max = 1- On the other hand, the fourth order differential equation 
( 13 . 1 p has the following oscillating solutions: 

S( v ) = C 1 (r j )e ik ^ dr ' + C 2 {r])e- lk J ~ dr] 



S 



ay/ 1 + 7r ay/l + f R 

We have used the WKB approximation under the condition % 2 <C k 2 . It is useful to define 

the effective growth factor f e g as f e s = dhi \ C\\/dN = d\n\C2\/N because the effective time 

evolution of the solution (13 .4p is determined by 6*1(77) or 62(77). In this case, the effective 

growth factor is given by 

which depends on the property of each model. If we consider some models of F(R) gravity, 
which does not satisfy the condition \H 2 fnR\ <C 1, the oscillating solutions can be growing 
ones. 



Case II 

Next, we consider the time evolution of the matter density perturbation in the case of 
(/c/K) 2 |/r| <C 1, I/rrI&Vo 2 S> 1 and % 2 <C k 2 . We can also obtain the following equation 



of matter density perturbation in this case by using Eq. (31) and the Appendix in 10j: 



5"" + H (3 + 0(H 2 /k 2 )) 5"' + U 2 (k 2 /H 2 + 0(1)) 5" 

4k 2 / W 



+ n 3 (k 2 /n 2 + o(i))5' -n 4 ^ 



jn 2 v h 2 

If we note to the terms proportional to k 2 , then we have 

d 2 6 (\ 3 \ d5 



i-^T? )+ 0(1)} 5 = 0. (3.6) 



dN 2 + {-2-2^)dN-^ 1+W ^ 6 -°- ™ 
In the matter dominant era, the equation 

d 2 5 1 d5 , . 

holds. The growth factors of this quasi-static solutions are / = — 1/4 ± a/33/4, which are 
same as those in Case I. While, the oscillating solutions are given as 

S( v ) = d-e lk ^ dv + C 2 -e- ikSd \ (3.9) 
a a 

where C\ and C 2 are arbitrary constants. The effective growth factor is expressed by 

/off = -I- (3.10) 

If \H 2 fnn\ <C 1, Eq. ( I3.10P is same as Eq. (j3.5]) . while if H 2 fnR is not negligible, it is different 
from that in Case I. 
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Case III 



If (k/H) 2 \fK\ <C 1, \fRR\k 2 /a 2 <C 1 and H 2 <ti k 2 , the fourth order differential equation 
of 5 is given by 

r/ W . / l2K 2 (-2 + H"/'H 3 )fRRR l-H'/H 2 nfaJ 2i 2\\ n, x >/> 

d + { + -2 + w/h> +0{u /x) r 5 f/ 

+X 2 1 (1 + 0(H 2 /x 2 )) S" + H(1 + 0(H 2 /x 2 )) S' + H 2 (^L - ^ + 0(H 2 /x 2 )) ^ = 0, 

(3.11) 



a 2 l-U'/U 2 

X =^Jh R 2-U»/W (3 ' 12) 
Noting to the terms proportional to x 2 i we obtain 

d 2 5 (\ 3 \ dS ( n H H . 



dN 2 \2 2 J dN \ H 2 H 3 

The solutions of Eq. (13. 13[) depend not only on /r and /rr but also on and /hrhr, 

because Friedmann-Lemaitre equations ( 12. 2j) and (I2.3P yield the following relation: 

2_^_ , 3r2 m fx fx 

# 2 # 3 ~ 2(1 + /*) + ^ 2 (1 + /«) 2W«(1 + /«) ' l ' J 

On the other hand, the oscillating solutions are given as: 

$(77) = c 1 e! f ' adN+i f xdn + C 2 ef f '* dN - i f xdn . (3.15) 
Here, Ci and C 2 are arbitrary constants, and the effective growth factor f cS is defined as: 

/cff = 1 " \W ln 1X1 - 2 ^N ln + 2-^3 - ( 3 - 16 ) 
We cannot further discuss the behavior of the solutions without assuming some models, 
because they also depend on higher derivatives of f(B). 

A. Model dependent analysis 

1. Hu and Sawicki's model 

We consider here the following form of the function f(R): 

m = -m* SI'" , (3.17) 
C2{R/m z ) n + 1 
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which was proposed by W. Hu and I. Sawicki 15]. If we use the best fit model[l6| to 
the Sloan Digital Sky Survey (SDSS) and the seven year data of the Wilkinson Microwave 
Anisotropy Probe (WMAP7), we can examine detailed time dependence of the matter den- 
sity perturbation. In this case, parameters of the model are given by n = 1.53, C\ = 10 3 ' 47 , 
and C2 = 10 2 28 , respectively, then this model is classified to Case III by assuming the scale, 
k ~ 0.1 — 0.01/iMpc" 1 . Therefore, the deviation from the ACDM model is a little at the 
perturbation order when we consider the quasi-static mode because R 3 /rrr and R 4 fRRRR 
are also negligible in these parameter sets. This result is consistent with the conclusions in 
The fast fluctuating mode will be considered a bit later. 



2. Starobinsky's model 



We next consider the Starobinsky's model 17] 



f(R) = XR 




(3.18) 



By the same procedure, we obtain the same result with Hu and Sawicki's model, i.e. the 
deviation from the ACDM model is so small that the observational verification is difficult 
when we consider the quasi-static mode. 

Such a result comes from the necessity to satisfy the local gravity tests. Namely, fitting 
the models of F(R) gravity to the general relativity in the solar system causes the deviation 
from the ACDM model a little also in the large scale. It is easy to understand where this 
result comes from by referring to Sec. IIB. That is to say, the condition \f-RR\k 2 /a 2 1 
induces this result. 



3. Exponential gravity model 



Next we consider the exponential type function of f{R) 18M20I]: 



f(R) = -cr(l-e 



-R/n 



(3.19) 



which has a possibility to deviate from the ACDM model, because of its sharp dependence 
of R. Some studies constrain the parameters in Eq. (13.191) 2ll |22J; however, there is no 
effective constraint on large c. If we consider c ~ 1000, we would find interesting behavior 
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of 5 of the quasi-static mode which comes from |-R 3 /Wr| ~ 1 , while \Rfn\, \R 2 /rr\ <C 1- 
However, it is difficult to cause the characteristic evolution of the matter density perturbation 
without changing the background evolution of the Universe, because \R 3 /rrr\ ~ 1 induces 
the changes in the background equations (12.21) and (I2.3p . 



4- The fast fluctuating mode 

We have seen that the F(R) gravity models which yield almost equivalent background 
evolution to that of the ACDM model cannot be also distinguished from the ACDM model 
at perturbation level when we consider the quasi-static mode of the matter density pertur- 
bation. So the other mode of the matter density perturbation, the fast fluctuating mode, is 
investigated in the following. The models which yield almost equivalent background evolu- 
tion to that of the ACDM model are classified as Case III, i.e. the case of (/c/"H) 2 |/rJ <C 1, 
|/rr|^ 2 / q2 1, and "H 2 <C k 2 . It is determined by the argument of \ whether or not the fast 
fluctuating solutions are oscillating solutions, because the behavior of the fast fluctuating 
solutions are understood from Eqs. (13.151) and (I3.16p . The argument of \ is found from the 
sign of fjiR because the relation, (1 — 'H'/'H 2 )/(2 — 'H"/'H 3 ) ~ 1, is held in this case. On the 



other hand, 
the scalaron 



;he condition > is imposed to the famous viable F(R) models to satisfy 
mass squared to be positive. Namely, x is re& l- Therefore, the fast fluctu- 



ating solutions always express oscillating behavior in those models. Next, it is important 
to investigate whether or not this oscillation is decaying oscillation. The expression of the 
effective growth factor (13.161) enables us to know the behavior of the oscillation. By using 
Eqs. fl22D and Q we have 



1 9 / H"\ n 2 f. 



2 2 V H 3 J a?f_ 



f eS c~ + ^[2-^-)^^. (3.20) 



' RR 

Here, "H"/"H 3 ~ 1/2 is held in the matter dominant era. Thus, it is determined by f^R and 
Jrrr that whether or not f e & is positive. If the values of parameters in 16( are applied to Hu 
and Sawicki's model and Starobinsky's model, we obtain / e ff < because of Jrr/ Jrrr < 0, 
i.e. it indicates these solutions are decaying oscillating solutions. This result does not 
strongly depend on the value of the parameters because the condition /rr/Jrrr < is 
derived from the property that Jr is approximately regarded as an inverse power law function 
of R in these models. At the same time, the oscillating solutions are also decaying solutions 
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in the exponential gravity model because /rr/ /rrr < is satisfied without depending on 
the values of parameters. 



The period of this oscillation is shorter than that of quintessence model ll|, as we can 
understand from the condition \f-RR\k 2 /a 2 <C 1. That is to say, this oscillation corresponds to 
the oscillation when the sound speed exceeds light speed in the fc-essence model. However, 
this oscillation does not mean superluminality in F(R) gravity. Therefore, there are not 
problems about causality. 

How to yield these solutions can be supposed from Eq. (I2.gp . The (ij), i ^ j element 
of the linearized Einstein equation simply yields $ = \1/ in the ACDM model and /c-essence 
model. However, there are solutions, satisfying $ = ^ or ^ ~ Jrr^" '/a 2 in F(R) gravity 
models. And the latter solutions correspond to the fast fluctuating mode. 

IV. CONCLUSIONS 

First, the time evolution of the matter density perturbation has been investigated by 
using approximations. When we use the sub-horizon approximation, it has been shown that 
the evolutionary behavior of the matter density perturbation depends on which condition 
is satisfied \fKR,\k 2 /a 2 ■€ 1 or |/rr|^ 2 /o 2 3> 1, as is well known. On the other hand, we 
may obtain the oscillating solutions by using the approximation d/drj ~ k instead of the 
sub-horizon approximation. And the variation of the matter density perturbation near the 
singularity has also been investigated as expanded applications of the approximate calcula- 
tions. 

Next, t he g rowth of the matter density perturbations has been investigated by using the 
results in [l0( and the case analysis. The first case is that the small scale approximation, 
H 2 -C k 2 , is more valid than others. The second case is that the approximation |/r| is small 
enough, and is more valid than small scale approximation, but |"H 2 /rr| is not very small. 
The third case is that approximations |/r| <C 1 and |H 2 /rr| <C 1 are more valid than the 
small scale approximation. In each of the cases, we have obtained characteristic solutions. 
Especially, the following behaviors have been shown: the quasi-static solutions in the first 
case grow faster than that of the ACDM model, and the fast fluctuating solutions in the 
third case have nontrivial dependence of the function f(R). The popular viable models 
correspond to the third case if the parameters in these models are determined to fit to the 
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observations of the background evolution of the Universe. And the quasi-static solutions 
are equivalent to that of the ACDM model. Moreover, the fast fluctuating solutions in 
those models are the decaying oscillating solutions. Thus, it is not possible to constrain 
these models from the matter density perturbation if the ACDM model indicates the true 
evolution of the background universe and the matter density perturbation. 

The conditions for the fast fluctuating mode to be decaying oscillating solutions are that 
of fun > and }'rrr < for the models , which yield almost equivalent background evolution 
to the ACDM model. So it can be said that imposing Jrr/Jrrr > enables us to make 
the model, which yields almost equivalent background evolution to the ACDM model, but 
yields different evolution of the matter density perturbation compared to the ACDM model. 
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